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Abstract 

It is known that the orientation of cellulose microfibrils within plant cell walls has an impor¬ 
tant impact on the morphogenesis of plant cells and tissues. Viewing the shape of a plant cell 
as a square prism or cylinder with the axis aligning with the primary direction of expansion 
and growth, the orientation of the microfibrils within the cell wall on the sides of the cell is 
known. However, not much is known about their orientation at the ends of the cell. Here 
we investigate the impact of the orientation of cellulose microfibrils within a plant cell wall 
at the ends of the cell by solving the equations of linear elasticity numerically. Three differ¬ 
ent scenarios for the orientation of the microfibrils are considered. The macroscopic elastic 
properties of the cell wall are obtained using homogenization theory from the microscopic 
description of the elastic properties of the cell wall microfibrils and wall matrix. It is found 
that the orientation of the microfibrils in the upper and lower parts of cell walls do not affect 
the expansion of the cell in the direction of its axis but do affect its expansion in the lateral 
directions. The arrangement of the microfibrils in the upper and lower parts of cell walls is 
especially important in the case of directed forces acting on plant cell walls and tissues. 

Keywords: biomechanics, plant modeling, homogenization, linear elasticity, plant cell wall 
microfibrils 

MSC subject classification: 35Q92, 74D05, 74Qxx, 92Bxx 


1 Introduction 

Knowing the influence of the microscopic molecular interactions and microscopic structure of 
plant tissues on the mechanical properties, development, and growth of plants is vital for the agri¬ 
culture and energy sectors. The mechanical properties of plant tissues are strongly determined 
by the mechanical properties of the cell walls surrounding plant cells and by the cross-linked 
pectin network of the middle lamella which joins individual cells together. Primary cell walls of 
plant cells, that are strong so as to resist a high internal hydrostatic pressure (turgor pressure) 
and flexible to permit growth, consist mainly of oriented cellulose microfibrils, pectin, hemicel- 
lulose, proteins, and water. The orientation, length, and high tensile strength of the microfibrils 
strongly influences the wall’s stiffness. Hemicelluloses form hydrogen bonds with the surface 
of cellulose microfibrils, which may effect the mechanical strength of the cell wall by creating 
a microfibril-hemicellulose network m- Pectin, once it is de-esterified and cross-linked with 
calcium ions, forms a gel within the primary cell wall and middle lamella and is hypothesized to 
be one of the main regulators of cell wall elasticity m 

Since the turgor pressure acts isotropically, it is the microstructure of the cell wall which 
determines the anisotropic expansion of plant cells. More specifically, it is the orientation of 
the cellulose microfibrils that influences the anisotropic expansion of the cell. Many plant cells, 
especially cells in plant roots and stem tissues, have a primary direction of expansion and less 
expansion takes place in the directions orthogonal to it, see e.g. mm- It is well-known that 
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cellulose microfibrils are parallel to the sides of primary cell walls and, particularly in young 
cells, perpendicular to the main direction of expansion [acZlEKIS]. For plant cells whose shape 
can be approximated by a prism or cylinder with the axis aligned with the primary direction of 
expansion, which is the case for root cells, this means that the microfibrils within the primary 
cell wall making up the sides of the cell are parallel to the sides and perpendicular to the axis 
of the cell. However, the orientation of the microfibrils in the cell wall at the ends of the cell 
does not appear to be known. Due to the important role that microfibrils play in the mechanical 
properties and expansion of the cell wall, knowing the orientation of the microfibrils everywhere 
is of vital importance. 

In this paper we investigate the affect of the orientation of the cellulose microfibrils in the 
upper and lower parts of cell walls on the deformation of the cell walls and plant tissues using 
numerical simulations. Modeling plant cells as a square prisms with rounded edges, we consider 
a part of a plant tissue represented by a central cell surrounded by cells on all sides. The 
primary cell wall and middle lamella are modeled as linearly elastic materials and on its internal 
boundary we specify a traction boundary condition corresponding to the turgor pressure. The 
cellulose microfibrils are arranged periodically within the cell wall, see e.g. m- The length scale 
of microfibrils (their diameter and distance between microfibrils) is much smaller than the scale 
associated with the thickness of the cell wall. This smaller length scale will be referred to as the 
microscale, while the scale associated with the dimensions of the cell wall is called the macroscale. 
To obtain the elastic properties of the primary cell wall we follow m and use homogenization 
theory to find an effective (macroscopic) elasticity tensor that depends on the orientation of the 
microfibrils on the microscale. It was observed experimentally that calcium-pectin cross-links 
influence mechanical properties of the cell wall matrix and of middle lamella, e.g. [2T]. The 
affect of the density of the calcium-pectin cross-links on the elastic properties of the cell wall 
are modeled through the Young’s modulus of the isotropic cell wall matrix. The microfibrils are 
assumed to be transversely isotropic. The effective elasticity tensors for cell walls are determined 
from the microscopic description of the mechanical properties of the microfibrils and cell wall 
matrix by solving numerically the corresponding ‘unit cell’ problems. Then using the macroscopic 
elasticity tensor for different microfibril orientations we solve numerically the equations of linear 
elasticity with different traction boundary conditions. The affect of the length of the cell in 
the direction of its axis on the deformation of the tissue is also investigated by considering two 
different cell lengths. 

We find that different configurations of orientations of microfibrils in the plane perpendicular 
to the main axis of cells has little effect on the expansion of the cells in the direction of its axis, 
however they do affect the expansion of the cell in the orthogonal directions. In general, we found 
that the expansion in the directions aligned with the microfibrils is less than the expansion in 
the directions orthogonal to the microfibrils. We also found that the expansion of the cell in 
the direction of its axis is smaller for shorter cells, which is in accord with Hooke’s law. If there 
are no applied forces and assuming the same turgor pressure in all cells, the expansion in every 
direction is smaller than in the case of additional forces acting on a plant tissue. The difference 
in the turgor pressure in the neighboring cells cause larger deformations in the directions parallel 
to the ends of cells, but the absolute value of the maximum of the deformation is negligibly 
affected by the orientation of the microfibrils in the upper and lower regions of the cell walls. 
The arrangement of the microfibrils in the upper and lower parts of the cell walls do have an 
impact on the elastic deformation of a plant tissue in the case where there are external forces or 
tissue tension in the directions parallel to the upper and lower parts of the cell walls. 

The outline of the paper is as follows. In Section we specify our model for plan tissue 
biomechanics. We consider the elastic deformation of the primary cell walls joined by middle 
lamella and the cell-inside is modelled by prescribing a turgor pressure. Next, in Section 
the results of our numerical simulations are presented. The numerical results are discussed in 
Section HI 


2 


2 Statement of the model 


In this section we present our model for the elastic deformations of a part of a plant tissue 
consisting of 27 cells connected by middle lamella. This section is divided into three parts: 
a description of the geometry of the domain, the presentation of the governing equations and 
boundary conditions, and the specification of the elasticity tensor on the domains representing 
the different parts of a plant cell wall and middle lamella. 

2.1 Geometry 

Our geometry is motivated by the structure of cells and tissues in young plant roots. We assume 
that the shape of a plant cell can be approximated by a square prism with rounded edges and 
consider a central cell surrounded by 26 cells, including diagonally adjacent cells. Choose a 
(xi,X 2 ,X 3 ) coordinate-system so that the origin is in the center of the central cell, the axes are 
parallel to the edges of the prism of the central cell, and the xs-axis is aligned with the axis of 
the central cell. Moreover, we consider 1 unit to be 1 /im. We consider the domain Q as depicted 
in Figurej^ the bounding box of which is (0,7.5) x (0,7.5) x (0,21.5). By reflecting f] over the 
planes xi = 0 , X 2 = 0 , and X 3 = 0 one obtains a domain that includes the central cell and parts 
of the 26 cells that surround it. The planes xi = 0, X 2 = 0, and X 3 = 0 will be called the planes 
of symmetry. 

A cross-section of 12 at a constant X 3 -value satisfying 0 < X 3 < 9.5 or 12 < X 3 < 21.5 is 
shown in Figure The regions with different orientations of the cellulose microfibrils and to the 
location of the middle lamella in this cross-section are shown in this figure. The cross-sections 
of 12 for a constant X 3 -value satisfying 9.5 < X 3 < 12 are not identical due to the rounded edges 
of the domain, see Figure [T] A cross-section of 12 at a constant X 2 -value satisfying 0 < X 2 < 2.5 
or 5 < X 2 < 7.5 is shown in Figure Once again, the domains with different orientations of the 
microfibrils and the location of the middle lamella are specified. Similar to the X 3 -direction, the 
cross-sections of 12 for a constant X 2 -value satisfying 2.5 < X 2 < 5 are not identical due to the 
rounded edges of the domain, see Figure A cross-section of the domain at a constant xi-value 
satisfying 0 < < 2.5 or 5 < < 7.5 is similar to Figure The thickness of 12 away from 

the junctions between two sections of the cell wall is 1.5 and the radius of all of the fillets is 0.5. 
The domain 12 is symmetric about the planes xi = 3.75, X 2 = 3.75, xs = 10.75, and xi = X 2 . 

The part of 12 contained within the box (0,7.5) x (0,7.5) X (10,11.5) is called the central 
region and represents the upper and lower parts of the cell walls. This region is divided into 
eight subdomains consisting of primary cell walls separated by a subdomain consisting of middle 
lamella. The eight subdomains consisting of primary cell walls are labeled 1-8 in Figure]^ The 
eight cells that make up 12 will be labeled according to which of these subregions they are in 
contact with. The eight subdomains consisting of cell walls are 3.6 /im in the xi and X 2 -directions 
and 0.6 /im in the X 3 -direction. The subdomain consisting of middle lamella separating these 
eight regions is 0.3 /im thick. To analyze the impact of the orientation of the microfibrils in the 
upper and lower regions of the plant cell walls on the elastic deformation of a plant tissue we 
will consider different microfibrils orientations within the subdomains 1 - 8 . 

2.2 Governing equations and boundary conditions 

The primary cell wall and the middle lamella are modeled as linearly elastic materials with 
different elastic properties. Let E be the elasticity tensor for the primary cell wall and middle 
lamella. The value of E = E(x) at any given point x G 12 depends on whether that point lies in 
the middle lamella or in the primary cell wall. Moreover, in the primary cell wall, the orientation 
of the cellulose microfibrils influences the elasticity tensor. This dependence will be specified in 
detail in the next subsection. 
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Figure 1: The domain ft consisting of parts of eight cells. By reflecting this domain over the 
planes xi = 0 , X 2 = 0 , and X 3 = 0 one obtains a domain that includes the central cell and parts 
of the 26 cells that surround it. 



2.5 fim 


0.5 fim 
0.6 /im 
0.3 /im 


Figure 2 : A cross-section of at a constant xs-value satisfying 0 < X 3 < 9.5 or 12 < X 3 < 21.5. 
All rounded corners in this flgure have a radius of 0.5 /im. The regions marked with 1 have 
cellulose microflbrils parallel to the xi-axis and regions marked with 2 have microflbrils parallel 
to the X 2 -axis. The region that is not marked is the middle lamella, which has no microflbrils. 
This cross-section is symmetric about the lines xi = 3.75, X 2 = 3.75, and xi = X 2 . 
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Figure 3: A cross-section of at a constant X 2 -value satisfying 0 < X 2 < 2.5 or 5 < X 2 < 7.5. 
All of the rounded corners have a radius of 0.5 /im. The regions marked with 2 have cellulose 
microfibrils parallel to the X 2 -axis, which is pointed into the page. The regions marked with c 
consist of the tops and bottoms of the cells and different microfibril orientations will be considered 
in these regions. The region that is not marked is the middle lamella, which has no microfibrils. 
This cross-section is symmetric about the lines xi = 3.75 and xs = 10.75. 
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Figure 4: A depiction of the subdomains within the central region of the domain. The eight 
regions labeled 1-8 are occupied by primary cell wall, and, hence, have cellulose microfibrils. 
Separating these regions is the middle lamella which is 0.3 /im thick. The regions with microfibrils 
are 3.6 /im on each side and are 0.6 /im thick. The fillets have a radius of 0.5 /im. The regions 
1-4 have larger X 3 -values than regions 5-8, and the xa-axis goes through the lower left corner of 
regions 3 and 7. 
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The boundary dQ of the domain can be split into the union of three sets: 


( 1 ) 

( 2 ) 

( 3 ) 


To = {x G dQ I = 0 or X 2 = 0 or X 3 = 0}, 

Tmax = {x ^ dQ \ xi = 7.5 or X 2 = 7.5 or xs = 21.5}, 

r/ = afi\(rour„,ax). 


The set T/ is the part of dQ in contact with the interior of the cells. A pressure boundary 
condition corresponding to the turgor pressure will be imposed on T/. On Tmax a tensile traction 
boundary condition will be specified. Finally, Tq is the part of the boundary of that lies on the 
planes xi = 0, X 2 = 0, or X 3 = 0 associated with the planes of symmetry. Thus, the displacement 
in the normal direction on Tq must be zero. 

Neglecting inertia and external body forces, the elasticity equation with these boundary 
conditions for the displacement u is 

div(Ee(u)) = 0 in fl, 

u • iz = 0 on To, 

< (Ee(u))i/ is parallel to u on Tq, (4) 

(Ee(u))l/ = fu on Tmax, 

3 Ee(u))i>'= —piv on F/, 

where e(u) = ^(Vu + Vu^) is the symmetric part of the gradient of the displacement and u is 
the exterior unit-normal to dQ. A unique solution of Q exists in [ 10 ] provided that 

/ G T^(rniax), P ^ and E satisfies the following conditions: 

1. |E| is bounded in 

2. There is a strictly positive a such that (a|Ap < A •E(x)A for all symmetric A G and 
X G fl. 

3. E possesses major and minor symmetries, i.e. 'Eijki = ^jikl — ^klij = ^ijlk- 


2.3 The elasticity tensor 

Next, we specify the elasticity tensor E on the domain Q. To do so, we must specify the elasticity 
tensor for the middle lamella and the primary cell wall for different microfibril configurations. 
The macroscopic elastic properties of the primary cell wall are derived from the microscopic 
description of the elastic properties of the cell wall matrix and microfibrils using homogenization 
theory. This requires the specification of the elastic properties of the cell wall matrix and the 
cellulose microfibrils. 

The cell wall matrix is isotropic 122 ], and so the elasticity tensor of the matrix 'Em is of the 
form 

E^vf A = 2/i7vf A + (tr A)l, 


where the Lame moduli fiM and Am are related to the Young’s modulus Em and Poisson’s ratio 
z/M through 


/iM(2/iM + 3 Am) 

XjM = -^- 


We take z/m = 0.3, which is common for biological materials, and Em = 5 MPa. This value is 
lower than the Young’s modulus measured for highly de-methylesterfied pectin gels considered 
in m since the pectin within the cell wall matrix is not fully de-esterfied. 

The cellulose microfibrils are not isotropic [3], so we assume that they are transversely 
isotropic and, hence, the elasticity tensor Ep for the microfibrils is determined by specifying 
five parameters: the Young’s modulus Ep associated with the directions lying perpendicular 
to the microfibril, the Poisson’s ratio i^pi characterizing the transverse reduction of the plane 


6 





perpendicular to the microfibril for stress lying in this plane, the ratio np between Ep and the 
Young’s modulus associated with the direction of the axis of the microfibril, the Poisson’s ra¬ 
tio up 2 governing the reduction in the plane perpendicular to the microfibril for stress in the 
direction of the microfibril, and the shear modulus Zp for planes parallel to the microfibriL A 
transversely isotropic elasticity tensor expressed in Voigt notation is given by 


/'q: 2 -h 0:5 

02 — as 

03 

0 

0 

A 

OL2 ~ Ct5 

02 + 05 

O3 

0 

0 

0 

as 

03 

OL\ 

0 

0 

0 

0 

0 

0 

04 

0 

0 

0 

0 

0 

0 

04 

0 

0 

0 

0 

0 

0 

05/ 


where for i = 1, 2, 3,4, 5, are related to the five parameters described above through 

Ep{l - ppi) 
npil - l^Fl) - 
Epup 

2nF{l - UFi) - 

EpUp2 

nF{l - UFi) - 

Zf, 

Ef 

2(1 + vfi) 

We assign these parameters the values 

Ep = 15,000 MPa, i^pi = 0.3, np = 0.068, i^P 2 — 0.06, Zp = 84,842 MPa, 

which are chosen based on experimental results [3] and to ensure that the elasticity tensor for 
the microfibrils is positive definite m 

We assume that the middle lamella is isotropic, with elasticity tensor and has a Young’s 

modulus of 15 MPa and Poisson’s ratio of 0.3. It is know from experiments that the density 
of calcium-pectin cross-links strongly influence the elastic properties of the cell wall matrix and 
middle lamella IZU- Thus, since in the middle lamella almost all pectin is de-esterified and the 
density of the pectin-calcium cross links is higher than in the cell wall matrix, where usually only 
70% of the pectin is de-esterified, we assume that the Young’s modulus for the middle lamella is 
three times larger than the Young’s modulus for the cell wall matrix. 

The cellulose microfibrils are arranged periodically within the cell wall matrix m and so 
standard techniques in homogenization theory, see e.g. m, yield a macroscopic elasticity tensor 
for a plant cell wall from the microscopic description of the mechanical properties of a cell wall 
on the level of a single mibcrofibril. In addition to the elastic properties of the microfibrils and 
cell wall matrix, the macroscopic elasticity tensor depends on the orientation of the cellulose 
mirofibrils. The components of this tensor are determined by solving unit cell problems, which 
have the form of the equations of linear elasticity and reflect the arrangement of the microfibrils. 

To specify the microstructure of a cell wall, consider the unit cell Y = (0,1)^ and let Ym 
and Yp represent the parts of Y occupied by the cell wall matrix and microfibrils, respectively, 
so that Ym and Yp are disjoint and Y = Ym ^Yp. Two configurations of microfibrils within Y 
are of primary interest. The first is when there is only one microfibril in Y occupying the set 

YF = {yeY\{y2- 0.5)^ + (ys - 0.5)2 ^ (5) 

and the other is when there are two microfibrils oriented in opposite directions and occupy 

YF = {yeY \iy 2 - 0.5)2 ^ ^ q_^ 252 or (yi - 0.5)2 ^ < 0.125^}, (6) 


ai = 

0^2 = 

0^3 = 
(a4 = 
0^5 = 
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Figure 5: A depiction of the unit cell Y with two configurations of microfibrils. Both unit cells 
are oriented so that the ?/ 3 -direction is pointing upward, (a) A picture of the unit cell with one 
microfibril occupying the set specified in (b) A picture of the unit cell with two microfibrils 
occupying the set specified in 


see Figure 1^ 

Then, the elasticity tensor Ey in Y is given by 


Ey(^) 


Em if 2/ e Tm, 
Ey if 2/ G ly, 


and can be extended T-periodically to all of M^. Consider a subdomain (7 of in which the 
cellulose microfibrils are arranged periodically with the orientation specified in T by ([^ or 
Let 6: be a small parameter associated with the distance between the cellulose microfibrils. The 
microfibrils of a plant cell wall are about 3 nm in diameter and are separated by a distance of 
about 6 nm, see e.g. la 01201, whereas the thickness of a plant cell wall is of the order of a 
few micrometers. To obtain the elasticity tensor for the part of the cell wall U with a periodic 
microstructure on the length scale of s defined by the structure of eY^ the periodic extension of 
Ey must be scaled appropriately. Namely, the elasticity tensor in U is given by 


E%x) = Ey for all x e U. 

Then homogenization theory yields a macroscopic elasticity tensor Ehom that describes a material 
whose behavior approximates the behavior of the cell wall with elasticity tensor E^ when 6: is 
very small uni. In our situation ps 10 Moreover, Ehom is given by 


^hom,ijkl i [^Y,ijkiiy) + ^Y,ijpqiy)ey{w^%q{y)]dy, 


( 7 ) 


where G iL^(y, M^) is the unique solution of 


j div(Ey(e(w^^) + b^^)) = 0 in T, 

1 fv ~ T-periodic, 

with 0 b^ + b^ 0 b^), where (b^, b^, b^) is the standard basis in M^. 

When Y is given by Q, the elasticity tensor given in Q will be denoted by as there are 
microfibrils in the xi and X 2 -directions, while when Y is given by the elasticity tensor given 
in 0 will be denoted by since the microfibrils are pointed in the xi-direction. Moreover, 

when ly is given by 0, then the microscopic elasticity tensor E^ depends only on two variables 






{x 2 ^xs) and the unit cell problem Q can be reduced to a two-dimensional problem [ 14 ]. To 
formulate this reduced problem, set Y = ( 0 , 1 )^ and 

Yf = ^Y\{y2- 0.5)2 ^ ^ q_252}, 

SO that y = ( 0 , 1 ) X y and Yp = ( 0 , 1 ) x Yjr. It can be shown that 


E- 


hom,ijkl ~ T. [^Y,ijkl{0^ v) + Ey^^jpg(0, ^)ey(w 
J Y 


ki\ 


Jpq 


{y)]dy, 


(9) 


with G being the unique solution of 

fdivy(Ey(0, y)(e(w^^) + b^^)) =0 


in Y, 


( 10 ) 


1 Iy ^y ~ y-periodic, 

where for a function w G H^{Y the differential operators and div^ are defined by 

/ 0 \ 

e(w) = Idy^^i dy^W 2 l{dy2^3 + dys^ 2 ) and div^w = dy^W 2 + dy^ws, 

Khdys^i U^y2^3 + dy3^2) J 

see e.g. M- Reducing the unit cell problem to two dimensions allows for the consideration of a 


higher resolution mesh when solving the problem (10) numerically. 

Besides considering the macroscopic elasticity tensor coming from microfibrils parallel to the 
xi-axis, we will also consider the macroscopic elasticity tensor generated by microfibrils that are 
arranged in other directions in the xiX 2 -plane. Given 9 G [—7r/2,7r/2], let denote the rotation 
about the X 3 -axis through the angle 0, so that 


cos 9 sin 9 0^ 
= ( — sin 0 cos 9 0 


0 


0 


1> 


The macroscopic elasticity tensor coming from a microstructure consisting of microfibrils 

aligned in the direction R^b^ is given by 

_ -p^ -dO -d6 -dO 11^1 

^hom,ijkl ~ ^ip^jq^kr^ls^hom,pqrs' 

So, for example, the macroscopic elasticity tensor coming from a microstructure with microfibrils 
parallel to the X 2 -axis is given by 

To summarize, the elasticity tensor E in the domain is different in different regions within 
the primary cell wall. In Figures and we specify the regions of the cell walls where the 
microfibrils are parallel to the xi-axis, i.e. E = regions of the primary cell wall 

where the microfibrils are parallel to the X 2 -axis, i.e. E = ' Within subregion i, for 

i — 1,...,8, of the central region, see Figure]^ the elasticity tensor E will be set equal to 
E^gj^, where different choices of E^^^^ associated with different microfibril configurations will be 
considered. Within the middle lamella, see Figures [2]-[^ there are no microfibrils and E = Eml- 
It follows from the properties of Em, E^d, and Eml that the macroscopic elasticity tensor E 
for the plant cell wall and middle lamella satisfies the conditions 1-3 mentioned at the end of 
Section 2.2 Hence the problem describing macroscopic elastic properties of plant cell walls 
connected by middle lamella is well-posed. 


3 Numerical results 

This section presents the results of the numerical simulations of the unit cell problems ^ and 
(10) necessary to calculate Ej^^^ and Ej^^^ and the simulations of the system @ for different 
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12.51 12.51 0 0 0 \ 

19.27 7.59 0 0 0 

7.59 19.27 0 0 0 

0 0 5.34 0 0 

0 0 0 9.30 0 

0 0 0 0 9.32/ 

Table 1: The effective elasticity tensor expressed in Voigt notation to two decimal places 

when the Young’s modulus of the matrix is 10 MPa. 


C^(IO) = 


/43333.24 
12.51 
12.50 
0 
0 
0 


configurations of cellulose microfibrils in the central region. All of the numerical simulations were 
done using FEniCS 13 El in]- This involved discretizing the domain using a nonuniform mesh 
and applying the continuous Galerkin method to solve the equations of linear elasticity. The 
resulting linear system was solved using the general minimal residual method with an algebraic 
multigrid preconditioner. 

3.1 Unit cell problems 

It was observed experimentally that the calcium-pectin chemistry influences the mechanical prop¬ 
erties of the cell wall matrix and middle lamella m Hence in general, the elastic properties of 
the cell wall matrix depend on the density of the calcium-pectin cross-links n and the microscopic 
elasticity tensor of the plant cell wall is a function of n. It was shown in m that under the 
assumption of an isotropic cell wall matrix, the macroscopic elasticity tensor Ehom corresponding 
to any microflbril conflguration is an affine function of the Young’s modulus of the cell wall 
matrix. From experiments 1221 , it is known that the Young’s modulus Em of the cell wall matrix 
is a function of the density of calcium-pectin cross-links n through the formula 

Em = 0.775n-h 8.08, (11) 

where Em has the units of MPa and n has the units of /iM. Thus, knowing the macroscopic 
elasticity tensor Ehom for two different values of Em determines the tensor for any value of Em- 
Then using 0 we obtain the macroscopic elasticity tensor for the cell wall for any calcium- 
pectin cross-links density n. This approach enables us to analyse the changes in the mechanical 
properties of plant cell walls and tissues in response to the dynamics of calcium-pectin chemistry 
and changes in calcium-pectin cross-link density, which will be the subject of future research. 

For the numerical simulations to obtain the macroscopic elasticity tensor we consider two 
values for the Young’s modulus: Em — 10 and Em — 20. Then using the fact that Ehom = 
IEhom(£’M) is an affine function we obtain Ehom for Em = 5. To And the unit cell Y was 

discretized by a mesh with 18,645,460 vertices with a higher density of vertices near the boundary 
between the cell wall matrix and the microflbrils. Using Voigt notation, the resulting effective 
elasticity tensor ^\^o^{Em) for Em = 10 and 20 are shown in Tables]^ andrespectively, to two 
decimal places. Using the symmetry of the microstructure it can be shown analytically that the 
entries of the matrices C^(IO) and C^(20) that are zero are exact and that some of the coefficients 
of the matrices C^(IO) and C^(20) are equal [15]. Speciflcally, for Em — 10 or 20, C^(£’m) 22 
and C^(£’m )33 should be equal, C^(£’m)i 2 and C^{Em)i 3 should be equal, and C^(£’m )55 and 
C^(£’m) 66 should be equal. The largest scale involved in the numerical computations of the 
macroscopic elasticity tensors is determined by the Young’s modulus of the microflbrils in the 
direction of the microflbrils and is equal to 220,000 MPa. Using this scale, the relative error 
associated with C^(£’m) 55 and C^(£’)66 not being equal is on the order of 10“^. 

Similarly, discretizing Y into a mesh with 11,750,289 vertices, the calculated effective elas¬ 
ticity tensor Ky^^{Em) for Em = 10 is shown in Tableand for Em = 20 is shown in Table 
to two decimal places, using Voigt notation. As before, it can be shown analytically that all 
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C^(20) 


/43352.40 24.07 
24.07 37.75 

24.07 14.89 

0 0 

0 0 

\ 0 0 


24.07 0 0 

14.89 0 0 

37.75 0 0 

0 10.44 0 

0 0 15.04 

0 0 0 


0 \ 

0 

0 

0 

0 

15.05/ 


Table 2: The effective elasticity tensor expressed in Voigt notation to two decimal places 

when the Young’s modulus of the matrix is 20 MPa. 



/10927.86 

99.60 

67.85 

0 

0 

0 \ 


99.60 

10927.69 

66.46 

0 

0 

0 


67.85 

66.46 

91.00 

0 

0 

0 

C^^(IO) = 

0 

0 

0 

186.83 

0 

0 


0 

0 

0 

0 

193.97 

0 


V 0 

0 

0 

0 

0 

352.58/ 


Table 3: The effective elasticity tensor expressed in Voigt notation to two decimal places 

when the Young’s modulus of the matrix is 10 MPa. 


of the entries of the matrices C^^(IO) and C^^(20) that are zero are exact and that some of 
the components of the matrices C^^(IO) and C^^(20) should be equal using the symmetry of 
the microstructure [15]. Specifically, for Em = 10 or 20, and C^^(£’m) 22 should be 

equal, C^^(£’m)i 3 and C^^(£^m )23 should be equal, and C^^(£’m) 44 and C^^(£^m) 55 should be 
equal. The largest relative difference between the components expected to be equal is of the 
order of 10“^. 

The results of this section allow us to compute the elasticity tensor for any Young’s modulus 
of the cell wall matrix, however in the following analysis we only consider the case where Em — 5 
MPa. 

3.2 Different boundary conditions and microfibril orientations in the central 
section 

Using the numerical results for the effective elasticity tensor for different microfibril orientations, 
in this section we consider different microfibril orientations in the eight subregions of the central 
section, see Figure]^ and different specifications of the turgor pressure p and tensile force / in 
problem @. 

We consider three different choices for p and /. To describe these, set po = 209 MPa, which 
is a common value for the turgor pressure [1]. 


C^^(20) 


A0943.35 

107.84 

75.25 

0 

107.84 

10943.18 

73.87 

0 

75.25 

73.87 

106.55 

0 

0 

0 

0 

191.43 

0 

0 

0 

0 

0 

0 

0 

0 


0 0 \ 

0 0 

0 0 

0 0 

198.56 0 

0 357.20/ 


Table 4: The effective elasticity tensor expressed in Voigt notation to two decimal places 

when the Young’s modulus of the matrix is 20 MPa. 
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Figure 6 : (a) A depiction of the orientation of the cellulose microfibrils in (Cl), (b) A depiction 
of the orientation of the cellulose microfibrils in (C 2 ). 


(BCl) Base case: p — Po and / = 5po . 

(BC2) No tensile tractions: p — Po and / = 0. 

(BC3) Different pressures, no tensile tractions: = P 4 = ps = ps = Po P2 — P3 — Pq — 

Py = 1.5po, where for i = 1,..., 8, is the pressure in cell i, and / = 0. 

For each of these boundary conditions we consider four different configurations of the microfibrils 
in the eight subregions of the center section. 


(Cl) In subregions 1 , 3, 5, and 7 the microfibrils are parallel to and in subregions 2 , 4, 

6 , and 8 the microfibrils are parallel to R“^/^b^.Thus, for i = 1, 3, 5, and 7, 
and for i = 2, 4, 6 , and 8 , Figure [^a). 

(C2) In subregions 2, 4, 5, and 7 the microfibrils are parallel to R^/^b^ and in subregions 1, 3, 
6 , and 8 the microfibrils are parallel to R“^/^b^. Thus, for i = 2, 4, 5, and 7, E^^^^ = ^ 

and for i = I, 3, 6 , and 8 , E^^^^ = See Figure [ gT). 

(C3) In all of the eight subregions the orientation of the microfibrils on the microscale are 
generated by the unit cell depicted in Figure [^b). Thus, E^^^^ = ^i^m. for i = I,..., 8 . 

(C4) There are no microfibrils in the center section. Instead, the central section consists of 
middle lamella and, hence, E^^^^ = E^l for i = I,..., 8 . 


The results of solving the system Q numerically using a mesh with 10,865,692 vertices with a 
higher density of vertices in the central region with the boundary conditions (BCI) and (BC 2 ) for 
the different configurations (CI)-(C4) are shown in TablesandFor the boundary condition 
(BC3), the system Q was solved on a mesh with 3,801,841 vertices with a higher density of 
vertices in the central region and the results are shown in Table A lower resolution mesh 
was used for (BC3) because the iterative solver failed to converge in 10,000 iterations when the 
higher resolution mesh was used. For each combination of boundary conditions and microfibril 
configurations the maximal displacement in the positive and negative xi, X 2 , and X 3 -directions 
are recorded to four significant figures. 

For configurations (Cl), (C 2 ), and (C4) and all boundary conditions the maximal displace¬ 
ments in the xi and X 2 -directions occur within the central region, while for configuration (C3) 
the maximal displacements in these directions occur on the sides of the cell wall. See Figure 
For the xs-direction, the maximal displacement occurs at the plane X 3 = 21.5. 
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(a) (b) 


Figure 7: A depiction of the displacements in the X 2 -direction for (BCl) with two different 
microfibril configurations: (a) configuration (Cl) and (b) configuration (C3) dark red regions 
correspond to the areas of maximal displacement in the positive X 2 -direction, while dark blue 
regions correspond to the areas of maximal displacement in the negative X 2 -direction. In these 
figures, the xi-axis is pointed out of the page and upward and the X 3 -axis is pointed to the right. 


(BCl) 

xi-direction 

X 2 -direction 

X 3 -direction 


negative 

positive 

negative 

positive 

negative 

positive 

(Cl) 

-0.05804 

0.1377 

-0.05691 

0.1380 

-3.080 X 10“^ 

2.010 

(C2) 

-0.05785 

0.2647 

-0.05696 

0.2659 

-8.263 X 10“^ 

2.009 

(C3) 

-0.05765 

0.05756 

-0.05639 

0.05790 

-3.779 X 10“^ 

1.965 

(C4) 

-0.05810 

0.2617 

-0.05682 

0.2617 

-4.496 X 10“^ 

1.981 


Table 5: The maximal displacement to four significant figures in the positive and negative xi, 
X 2 , and X 3 -directions for boundary condition (BCl) and microfibril configurations (C1)-(C4). 


3.3 Smaller cells 


Besides considering the situations mentioned in the previous subsection, we also consider the 
case where the cells are smaller. Namely, we consider cells 10 /im smaller than those described in 
Section 2J_ so that the length of the domain in the X 3 -direction is 11.5 /im. Looking at Figure 
this means that the 10 /im measurements are decreased to 5 /im. Moreover, the boundaries 0 
and 0 must be replaced with 


Fo = {x G dQ I = 0 or X 2 = 0 or X 3 = 5}, (12) 

Tmax = {x G I = 7.5 or X 2 = 7.5 or xs = 16.5}. (13) 

For the case of smaller cells, which will be referred to as (SM), we consider boundary condition 
(BCl) and configurations (C1)-(C4). The results of solving the system 0 numerically using a 
mesh with 9,490,948 vertices with a higher density of vertices in the central region are shown in 
Table [H 


4 Discussion 

The data in Tables [5}|^ tell us several things about the affect of the presence and orientation of 
the cellulose microfibrils in the central region, i.e. the upper and lower ends of cell walls. First 
of all, they have little affect on the expansion of the cells in the X 3 -direction, as can be seen 
from looking at the last two columns in these tables. The cell wall is able to expand more in the 
directions perpendicular to the directions of the microfibrils since the microfibrils are much stiffer 
than the cell wall matrix and middle lamella. Thus, changing the microfibril orientation within 
the xiX 2 -plane has little affect on the displacement in the X 3 -direction. However, the expansion 
in the xi and X 2 -directions are affected. In particular, for (BCl) and (SM) when the microfibrils 
are arranged in the configuration (C3), the displacement in the positive xi and X 2 -directions 
are 2/5 of those for configuration (Cl) and 1/5 of those for configurations (C2) and (C4). In 
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(BC2) 

xi-direction 

X2-direction 

X3-direction 


negative 

positive 

negative 

positive 

negative 

positive 

(Cl) 

-0.02037 

0.01509 

-0.01968 

0.01491 

-6.208 X 10“^^ 

0.1183 

(C2) 

-0.02034 

0.01451 

-0.01969 

0.01478 

-1.881 X 10“^*^ 

0.1184 

(C3) 

-0.02035 

0.01425 

-0.01967 

0.01415 

-6.416 X 10“^*^ 

0.1209 

(C4) 

-0.02031 

0.01640 

-0.01967 

0.01628 

-1.036 X 10-1“ 

0.1212 


Table 6: The maximal displacement to four significant figures in the positive and negative xi, 
X 2 , and X 3 -directions for boundary condition (BC2) and microfibril configurations (C1)-(C4). 


(BC3) 

xi-direction 

X2-direction 

X3-direction 


negative 

positive 

negative 

positive 

negative 

positive 

(Cl) 

-0.005638 

0.2041 

-0.05627 

0.03504 

-2.405 X 10-1“ 

0.1449 

(C2) 

-0.007259 

0.2017 

-0.05551 

0.03476 

-1.299 X 10-1“ 

0.1451 

(C3) 

-0.007191 

0.1990 

-0.05243 

0.03441 

-5.5574 X 10-11 

0.1488 

(C4) 

-0.005591 

0.2045 

-0.05604 

0.03508 

-2.789 X 10-1“ 

0.1477 


Table 7: The maximal displacement to four significant figures in the positive and negative xi, 
X 2 , and X 3 -directions for boundary condition (BC3) and microfibril configurations (C1)-(C4). 


configuration (C3) there are microfibrils oriented in both the xi and X 2 -directions within the 
central region and it is expected that for this configuration there would be less expansion in both 
directions. The difference in the maximal deformations for different microfibril configurations 
in the central region is less noticeable for the (BC2) and (BC3) boundary conditions. Hence 
these results indicate that the orientation of microfibrils in the central section has an important 
impact on the deformation of plant cell walls and tissues in the case of tensile traction boundary 
conditions. 

Comparing Tables and we see that the presence of the tensile traction boundary condition 
causes the displacements in the positive directions to increase by an order of magnitude. This is 
not surprising as the presence of more forces causes larger displacements. These tables also show 
that in the absence of a tensile traction boundary condition, the case (BC2), the displacement 
in the negative xi and X 2 -directions is greater than the displacement in the positive xi and 
X 2 -directions. This is not the case for (BCl) and (SM) because of the applied tensile forces 
in the positive directions. For (BC2), the reason that the absolute value of the displacement 
in the negative directions are greater than the maximal displacements in the positive xi and 
X 2 -directions is the result of several factors. Let us focus our discussion on the xi-direction 
as the X 2 -direction is similar. The first thing we notice is that there are microfibrils in the 
xi-direction, so any displacement in this direction is greatly hindered. Next, we consider the 
boundary conditions on T/ and T^ax and the fact that the turgor pressure p is balanced with 
respect to the positive and negative directions. The surfaces where the pressure is pointing in 
the positive xi-direction are closer to a surface in which a no displacement in the xi-direction 
boundary condition is imposed than the surfaces where the pressure is pointing in the negative 
xi-direction. Thus, the no displacement in the xi-direction hinders the displacement in the 
positive xi-direction more than the displacement in the negative xi-direction. This reasoning 
cannot be applied to the X 3 -direction because there are no microfibrils pointing in X 3 -direction. 

Comparing Table and Table we can see the effect of increasing the pressures in some of 
the cells. First, notice that the displacements in the positive directions are larger in the case 
(BC3) than in the case (BC2). This is because in (BC3) the pressure in cells 2, 3, 6, and 7 
is greater than in the (BC2) case. Also notice that in the case (BC3) the displacement in the 
positive xi-direction is greater than the displacement in the positive X 2 and X 3 -directions, which 
is caused by the the position of the cells with the larger pressure. Namely, there is a pressure 
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(SM) 

xi-direction 

X 2 -direction 

X 3 -direction 


negative 

positive 

negative 

positive 

negative 

positive 

(Cl) 

-0.06053 

0.1377 

-0.06043 

0.1379 

-3.583 X 10-*^ 

1.066 

(C2) 

-0.05811 

0.2665 

-0.05788 

0.2657 

-3.942 X 10“^ 

1.058 

(C3) 

-0.05744 

0.05842 

-0.05692 

0.0559 

-7.333 X lO-*^ 

1.011 

(C4) 

-0.06285 

0.2619 

-0.06275 

0.2621 

-2.214 X 10"^ 

1.041 


Table 8: The maximal displacement to four significant figures in the positive and negative xi, 
X 2 , and X 3 -directions for the case of smaller cells for microfibril configurations (C1)-(C4). 


difference between the cells that are aligned in the xi-direction. The impact of the presence and 
orientation of microfibrils on the displacement is relatively small. 

Finally, looking at Tables and one sees that the only significant difference in the case 
of small cells is in the displacement in the xs-direction—there is twice the displacement for the 
larger cells than for the smaller cells. This is in accord with Hooke’s law, which tells us that the 
elongation of an elastic bar under an applied load is a linear function of the length of the bar. 

To conclude, the orientation of the microfibrils in the upper and lower parts of plant cell 
walls have no effect on the elongation of the cells, but will influence their radial expansion and 
growth. The results in Table show that different pressures in neighbouring cells, which can be 
observed during the growth process, influence the direction of the maximal displacement (here 
the maximal displacement in the xi-direction is due to pressure distributions). It follows from 
our results that only in the case of directed tensile forces applied to plant cells and tissues will 
the orientation of the microfibrils in the lower and upper parts of cell walls play a role. Hence, 
for cells where the main acting forces are turgor pressure, the orientation of the microfibrils in 
the lower and upper parts of the cell walls is not essential and cells will choose the most energy 
efficient way to orient the microfibrils in these parts of the cell walls. However in the parts of the 
tissues where there is a strong directed tissue tension, the importance of the orientation of the 
microfibrils in the upper and lower parts of the cell walls may be important. In our studies we 
assumed that the microfibrils on the sides of the cell wall are arranged in fixed rings around the 
cells without considering possible sliding of the microfibrils during the expansion. The affect of 
the sliding of the microfibrils on the deformation of plant cells and tissues in combination with 
different arrangements of microfibrils in the upper and lower parts of the cell walls will be the 
subject of future studies. 
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